NMR multiple quantum coherences in quasi-one-dimensional spin systems: 
Comparison with ideal spin-chain dynamics 
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The spins in a crystal ol fluorapatite have often been used to experimentally approximate a 
one-dimensional spin system. Under suitable multi-pulse control, the nuclear spin dynamics may 
be modeled to first approximation by a double-quantum one-dimensional Hamiltonian, which is 
analytically solvable for nearest-neighbor couplings. Here, we use solid-state nuclear magnetic res- 
onance techniques to investigate the multiple quantum coherence dynamics of fluorapatite, with an 
emphasis on understanding the region of validity for such a simplifled picture. Using experimental, 
numerical, and analytical methods, we explore the effects of long-range intra-chain couplings, cross- 
chain couplings, as well as couplings to a spin environment, all of which tend to damp the oscillations 
of the multiple quantum coherence signal at sufflciently long times. Our analysis characterizes the 
extent to which fluorapatite can faithfully simulate a one-dimensional quantum wire. 

PACS numbers: 03.67.Hk, 03.67.Lx, 75.10.Pq, 76.90.-|-d 



I. INTRODUCTION 

Low-dimensional quantum spin systems are the sub- 
ject of intense theoretical and experimental investigation. 
From a condensed matter perspective, not only do these 
systems provide a natural setting for deepening the ex- 
ploration of many-body quantum coherence properties 
as demanded by emerging developments in spintronics 
and nanodevices [J, d, 0] , but the ground states of one- 
dimensional (ID) conductors provide insight into the so- 
lution of the one-band Hubbard Hamiltonian Q . From 
a quantum information perspective quantum spin 
chains have been proposed as quantum wires for short- 
distance quantum communication, their internal dynam- 
ics providing the mechanism to coherently transfer quan- 
tum information from one region of a quantum computer 
to another (see also Q for a recent overview) . Perfect 
state transfer, in particular, has been shown to be the- 
oretically possible by carefully engineering the couplings 
of the underlying spin Hamiltonian. A number of efforts 
are underway to devise protocols able to achieve reliable 
quantum information transfer under more realistic condi- 
tions - bypassing, for instance, the need for initialization 
in a known pure state Q, explicitly incorporating the 
effect of long-range couplings [ o l, [lOl . [ll|, or exploiting 
access to external end gates [12|, • Still, few (if any) 
physical systems can meet the required constraints, and it 
is likely that quantum simulators will be needed to exper- 
imentally implement these schemes. Of course, quantum 



simulators will in turn allow us to probe a much broader 
range of questions encompassing both quantum informa- 
tion and condensed matter physics [ij]- Optical lattices 
have shown much promise in simulating quantum spin 
systems [l^. Among sohd-state devices, coupled spins 
in apatites have recently enabled experimental studies of 
ID transport and decoherence dynamics [igI. ITtI. [TsI. [4^ . 

Fluorapatite (FAp) has long been used as a quasi-lD 
system of nuclear spins. Lowe and co-workers character- 
ized the nuclear magnetic resonance (NMR) line shape of 
FAp [l^, [2^ , and described the dipolar dynamics of the 
free induction decay in terms of the ID XY model [2l|. 
Cho and Yesinowski investigated the many-body dy- 
namics of FAp under an effective double-quantum (DQ) 
Hamiltonian, and showed that the growth of high-order 
quantum coherences was distinctly different from that 
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obtained in dense 3D crystals 
cal standpoint, FAp provides a rich testbed to explore the 
controlled time evolution of a many-body quantum spin 
system. The DQ Hamiltonian is analytically solvable 
in the tight-binding limit, where only nearest neighbor 
(NN) couplings are present [3, [2^, [2^ . Previous work 
showed that the implementation of a DQ Hamiltonian in 
the FAp system using coherent averaging techniques is 
a promising tool for the study of transport in quantum 
spin chains. We demonstrated, in particular, that the 
DQ Hamiltonian is related to the XY-Heisenberg Hamil- 
tonian by a similarity transformation, and that it is pos- 
sible to transfer polarization from one end of the chain 
to the other under the DQ Hamiltonian [13]. In fact, 
the signature of this transport shows up in the collective 
multiple quantum coherence (MQC) intensity of the spin 
chain. Experimentally, it has also been shown that it is 
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possible to prepare the spin system in an initial state in 
which the polarization is localized at the ends of the spin 
chain [l6l |. pavin g th e way towards achieving universal 
quantum control [27| . 

Since the mapping between the experimental system 
and the idealized model 0, [i3| is not perfect, an es- 
sential step forward is to address where and how this 
model breaks down, which constitutes the main aim of 
this paper. In particular, we systematically examine the 
viability of using NMR investigations of FAp as a test- 
bed for ID transport, by relying on a combination of 
experimental and numerical methods. We first examine 
the effects on the relevant observables of experimental 
errors introduced during the implementation of the DQ 
Hamiltonian, which arise due to higher-order terms in the 
average Hamiltonian describing the effective spin evolu- 
tion. We also examine errors introduced in some state 
initialization sequences due to the restriction of the con- 
trol fields to collective rotations. Since the FAp crystal is 
in reality a three-dimensional (3D) lattice, we next inves- 
tigate in detail how the spin dynamics is affected by the 
presence of longer-range couplings, both within a single 
chain and between adjacent spin chains. 

The content of the paper is organized as follows. We 
describe the quasi-lD spin system of FAp in Sec. [Tll in- 
cluding the evolution in the absence of control as well 
as the dynamics under suitable pulse sequences. In the 
same section, we also discuss the system initialization 
and the readout of the experimental MQC signal. Sec- 
tions mil and IIVI present both experimental and numer- 
ical results of MQC dynamics, and are the core of the 
paper. By comparing the numerical results with the an- 
alytical predictions available in the limiting case of a DQ 
Hamiltonian with NN couplings, we evaluate the effect 
of high-order average Hamiltonian terms, next-nearest- 
neighbor (NNN) couplings, and cross-chain couplings be- 
tween multiple chains. Our findings are summarized in 
Sec. |Vl Appendix |X] presents technical background on 
the relevant numerical methodology, whereas we also in- 
clude in Appendix [B] a description of finite size effects as 
found in simulations, and in Appendix [Cl a discussion of 
an alternative chaotic model for the spin bath. 



II. PHYSICAL SYSTEM AND EXPERIMENTAL 
SETTINGS 

A. Spin Hamiltonian of fluorapatite 



tization field, leading to a Hamiltonian of the form: 



N 
j<£ 



(1) 



Here, cr" (a — x, y, z) denotes the Pauli matrices of the 
jth spin and bj^ = (/^o/167r)(7^/i/r|^)(l — 3 cos^ Oj^), with 
/io the standard magnetic constant, 7 the gyromagnetic 
ratio of fluorine, Tji the distance between nucleus j and 
£, and Oji the angle between rj^ and the z-axis. The ge- 
ometry of the spin system is reflected in the distribution 
of the bji couplings. 

The FAp crystal has a hexagonal geometry with space 
group PGs/m (see Fig. [T]a). The dimensions of the 
unit ceU are D = 9.367 A and c = 6.884 A. The ^^F nuclei 
form linear chains along the c-axis, each one surrounded 
by six other chains. The distance between two intra-chain 
^^F nuclei is d = c/2 = 3.442 A and the distance between 
two cross-chain ^^F nuclei is D. The largest ratio between 
the strongest intra- and cross- chain couplings (w 40) is 
obtained when the crystalline c-axis is oriented parallel to 
the external fleld. Thus, to a first approximation, in this 
crystal orientation the 3D ^^F system may be treated as a 
collection of many identical ID chains. For a single chain 
oriented along z, we have bji = — (^o/''r)(7^ft/c'^|j — £p). 
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FIG. 1: a) Unit cell of the fluorapatite crystal [Ca5(P04)3F], 
highlighting the geometry of the fluorine chains (in red), b) 
NMR pulse sequence for the creation and detection of MQC. 

In reality, naturally occurring defects in the sample 
(such as vacancies or substitutions [l^]) cause the chains 
to be broken into many shorter chains. Here we model 
the system as an ensemble of (approximately) indepen- 
dent and equivalent chains with finite length. Such a 
simplified description is necessary to obtain a computa- 
tionally tractable model. 



We consider a single crystal of FAp [Ca5(P04)3F] at 
room temperature, placed in a strong external magnetic 
field along the z-direction that provides the quantiza- 
tion axis for the nuclear spins. It is possible to truncate 
the magnetic dipolar interaction among the spins in this 
strong field, keeping only the secular terms. The result- 
ing secular dipolar interaction [2^ among N ^^F nuclear 
spin- 1/2 is anisotropic due to the presence of the quan- 



B. Control capabilities and effective dynamics 

1. Unitary control 

Unitary control is obtained by applying (near)resonant 
radio-frequency (rf) pulses to the spin system. FAp con- 
tains ^^F and '^-^P spins-1/2, both of which are 100% 
abundant. Moreover, in an ideal crystal, all the ^^F spins 
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are chemically equivalent, as are all the ^^P spins. As a 
consequence, all rf control pulses are applied collectively 
to all the spins. 

In NMR, the term MQC refers to coherences between 
two or more spins. When the system is quantized along 
the z-axis, a quantum coherence of order n is asso- 
ciated to the transition between two states \mi) and 
1 771,2), such that the difference of the magnetic moment 
along z of these states (toi — 771,2) oc n. That is. mul- 
tiple quantum coherences of order n describe states like 
|'Ti2)("^i|i or elements in the density matrix that cor- 
respond to a transition between these two states js^ . 
Quantum coherences can also be classified based on 
their response to a rotation around the z (quantization) 
axis. A state of coherence order n acquires a phase pro- 
portional to n under a z-rotation. Multiple quantum 
NMR techniques [H, [H, [H, [H, [s^ have enabled re- 
searchers to probe multi-spin processes, and gain insight 
into the many- body spin dynamics of dipolar-coupled 

solids m, m, so, m, H S, El ■ 

To study the MQC dynamics of the spin system, we 
typically let it evolve under the DQ Hamiltonian 

= (2) 

with af = (crj ± icrj')/2. Following [H,!!!], we utilize a 
16-pulsc cycle applied on-rcsonance with the ^^F Larmor 
frequency to implement the DQ Hamiltonian to lowest 
order in AHT description [H, |4^, [4^ . A key feature of 
this sequence is that the fluorine-phosphorus dipolar in- 
teraction is decoupled, which makes it possible to ignore 
the presence of the "^^P spins in the rest of this paper. 
The dynamics during the pulse sequence can be written 
in terms of an effective Hamiltonian TLdq^ 

= e±'^°«*, (3) 

where T denotes time-ordering operator, h = I, and 
H^''^\t) is the time-dependent Hamiltonian describing 
the rf-pulses along the x- (or y-) axis (whereby the cor- 
responding ± sign in front of the effective Hamiltonian) . 

2. Initialization capabilities 

The spin dynamics under the DQ Hamiltonian depends 
critically on the initial state in which the system is pre- 
pared. Here, we focus our attention on two choices of di- 
rect experimental relevance [l6| . One is the equilibrium 
Zeeman thermal state, which is obtained at the thermal 
equilibrium in a strong external magnetic field (Bq = 7 T 



in our experiments) at room temperature. The thermal 
state can be expressed as 

PthiO) oc exp(-ecr^) « 1 - ea^, (4) 

where cr^ = J2j '^j s-^^d e = jBo/kBT, with fc^ the Boltz- 
mann constant and T the temperature (e « 10~^ at room 
temperature for FAp). In line with standard NMR prac- 
tice, we consider only the evolution due to the compo- 
nent proportional to e, pth(O) = cr^, since the identity 
matrix does not evolve or contribute to the MQC signal 
under the assumption of unital dynamics. The second 
initial state that is experimentally available is a mixture 
of states where only a spin at the extremities of the chain 
is polarized, which can be formally represented as 

Pe,ld{0) = C^I+CTn: (5) 

where spin 1 and spin TV are located at the two ends of 
the spin chain p^ . We refer to this as the end-polarized 
state. A description of the method used to create this 
state is given in Sec. IIV A2I 

3. Readout capabilities 

In an inductively detected NMR experiment (in which 
a coil is used to measure the average magnetization) , the 
observed signal is S{t) = ((a~{t)) = (^Tr {(T^p(i)}, where 
(7^ = aj and C is a proportionality constant. The 
only terms in pit) that yield a non-zero trace, and there- 
fore contribute to S{t), are angular momentum opera- 
tors such as CT^, which are single-spin, single-quantum 
coherences. Thus, in order to characterize multi-spin dy- 
namics, it is necessary to indirectly encode the signature 
of the dynamics into the above signal. This is precisely 
what is done in standard NMR MQ spectroscopy, using 
an evolution-reversal experiment [321 • The density oper- 
ator at the end of an MQ experiment is given by 

Pf U\jqU^^o\UmqPiUIiqUI^^^Umq, (6) 

where Umq = exp(— iTiugf;), and [/ovoi determines the 
nature of the information encoded (see Fig. [TJb). 

In our experiment, we are interested in the evolution of 
MQC under the DQ Hamiltonian, thus wc measure the 
signal as we systematically increase t. In order to encode 
information about the distribution of the MQC, we ap- 
ply a collective rotation about the z quantization axis, 
Uevoi = exp(— z0(Tz/2). Then, to extract the coherence 
order distribution, the measurement is repeated while 
incrementing from to 27r, in steps of 5(j) = l-n jlK, 
where K is the highest order of MQC encoded. The sig- 
nal acquired in the fc-th measurement is then S^(t) = 
Tr[/9'"'(i)(7z], where p^it) is the density matrix evolved 
under the propagator Uk{t) = e>''><-oQi(>-ikS<l,a,ji^-mDQt ^ 
and we have assumed that is the experimental observ- 
able. In practice, we use either ai: /2 pulse or a solid echo 
[47t to read out the signal at the end of the experiment. 
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Fourier-transforming the output with respect to 
the coherence order intensity: 



Jn{t) 



K 

E 



' yields 



(7) 
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Note that since the initial states we consider are popula- 
tion terms in the z-basis, the final states at the end of the 
evolution-reversal experiment are also population terms 
(hence our use of the observable cr^). 



III. MULTIPLE QUANTUM DYNAMICS: 
SIMPLE MODEL AND EXPERIMENT 

A. Ideal spin-chain dynamics 

The fact that the evolution of the ID spin chain under 
a DQ Hamiltonian is exactly solvable in the tight binding 
limit [I^, [2^, [2^ provides a useful starting point for the- 
oretical analysis. Hereafter, we shall refer to this model 
as the analytical model. Moreover, the DQ Hamiltonian 
is related to the XY Hamiltonian by a similarity trans- 
formation that inverts every alternate spin in the chain 
[itI . ■ Besides using the analytical results to calibrate 
our numerical methods (see Appendix |A| , we will also 
investigate the effect of long-range interactions beyond 
the NN limit, by comparing numerical and analytical re- 
sults. For convenience, we set the NN coupling strength 
in the DQ Hamiltonian of Eq. ([2]) 6 = fci2 = 1, so that 
time shall be measured in units of 1/6 henceforth (unless 
explicitly stated otherwise). 

For both the thermal and the end-polarized initial 
state, only zero and DQ coherences are predicted by 
the analytical model. Specifically, for the thermal ini- 
tial state, the normalized intensities are 



k 

= ;-^5]sin2(46tcos^fc). 



(8) 



where as before N is the number of spins in the chain and 
ipk = k7T/{N + 1). For the end-polarized initial state. 



Jo'" it) = 

In both Eqs. 
that J0 + 2J2 



N 
TV 



— sin^ ( V-fc ) cos^ {ibt cos V-fc ) , 
fc 

^^sin2(V'fc)sin2(46icosV'fc)- (9) 



the normalization is chosen such 



B. Experimental results 

The experiments were performed in a 7 T magnetic 
field using a Bruker Avance Spectrometer equipped with 
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FIG. 2: (Color online) Experimental data (OQ circles, 2Q- 
stars, 4Q squares), analytical results (red dash-dotted lines), 
and least square fit (blue dashed lines) of experimental data 
to numerical calculations under the DQ Hamiltonian with NN 
and NNN couplings of strength b and b/8, respectively. Error 
bars for the experimental data were estimated from the stan- 
dard deviation of the odd quantum coherences from zero. Nu- 
merically, two free parameters, the time origin to and the time 
scale 1/b, are adjusted to obtain the best fit to the experimen- 
tal data for Jq. Top: Thermal initial state with chain length 
A*' — 18. Fitting parameters to — 9.6/js and 1/6 = 101.8/is. 
Bottom: End-polarized initial state wth chain length TV = 19. 
Fitting parameters to — 0.2/is and 1/6 = 100.8/is. 



a home-built probe. The ^^F frequency is 282.37 MHz. 
Experimentally measured MQC data are shown in Fig. [21 
along with analytical predictions and simulation results 
under the DQ Hamiltonian with NN and NNN couplings. 
Both the time origin to and the coupling strength b were 
used as fitting parameters in order to minimize the square 
of the difference between the experimental and numerical 

data, that is, J2r iJo^'^^iU) - JS^'^'lHU - ta)]f- For the 
thermal initial state data (upper row in Fig. [2]), the 7r/2 
pulse length was 1.05 fis^ the inter-pulse delay A was 
varied from 2.9 to 5 fis, and the number of loops was 
increased from 1 to 7. We set K = 12, and incremented 
the phase in steps of 27r/24 to encode the MQCs. For the 
end-chain initial state data (lower row in Fig. [2]), the 7r/2 
pulse length was 0.93 ^s, the end-state preparation time 
ti = 30.3 fis, the inter-pulse delay A was varied from 2.9 
to 7.3 /iS, and the number of loops increased from 1 to 
8. We set K = 16, and incremented the phase in steps 
of 27r/32 to encode the MQCs. In both cases, the recycle 
delay was 300s and a solid echo sequence with an 8-step 
phase cycle was used to read out the signal intensities at 
the end of the experiment. 

The experimental data are normalized at every time 
step, such that Jq -I- 2(J2 -I- J4) = 1 (using the fact 
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that J_„ = J„). The intensities of the odd MQCs (not 
shown in Fig. [5]) turn out to be negligibly small. At 
short times (less than sa 0.2 ms), Fig. [2] indicates that 
fourth- and higher-order even-MQCs are also negligible. 
However, the four-quantum coherence signal contributes 
significantly at longer times. In 3D systems, including 
both plastic crystals such as adamantane [4a| and rigid 
crystals such as the cubic lattice of -'^^F spin in CaF2 
[iol . l43j , very high coherence orders are seen to develop 
over a time scale less than a millisecond, with no appar- 
ent restriction on the highest order reached. In contrast, 
the fact that the MQ intensities are restricted to the zero- 
and DQ-coherences, and that the higher-order terms only 
grow relatively slowly during the whole time domain we 
explored, are strong indications of the ID character of 
the spin system. At the same time, the appreciable in- 
tensity of the four-quantum coherence at long evolution 
times clearly indicates that the analytical model (which 
predicts only zero- and DQ- coherences for both the ther- 
mal and end-polarized initial states) becomes inadequate 
to accurately describe the real system. 

Note that in the simulations, the maximum compu- 
tationally accessible chain length was = 21 spins. 
Though the fits included in Fig. 2 use 18-19 spins, it 
is important to realize that sensitivity of the dynamics 
to the precise value of N develops only at sufficiently long 
times (as the effect of the finite chain boundaries mani- 
fest - see Appendix [B|, where the accuracy of the simple 
model used to make the estimate becomes itself limited. 



A. Errors due to limited control 

1. High-order terms in Average Hamiltonian Theory 

As mentioned, experimentally the DQ Hamiltonian ^ 
is obtained as the zeroth order average Hamiltonian of a 
multiple-pulse sequence. Since the 16-pulse cycle used 
in the experiment is time-symmetric [42|, all odd-order 
corrections are zero, and the leading error term is of 
second-order in HTY'^'TcH , both when considering ideal 
and finite-width pulses. The contributions to the effective 
Hamiltonian from higher-order terms may be estimated 
by comparing the single-cycle MQC signal computed us- 
ing the exact DQ Hamiltonian, and using the dipolar 
Hamiltonian ^ interspersed with rf pulses, respectively. 
Assuming ideal instantaneous pulses, we verified numer- 
ically that for the system of interest such contributions 
arc small provided that the cycle time Tc < 4 (sec Fig. |31 
Inset). In the experiment, we thus employed multi-cycle 
sequences, in order to extend the region of validity of the 
DQ model. 

In order to determine how well we implemented the 
evolution reversal experiment described in Sec. II, we 
performed a series of experiments that measured the 
overlap between the initial and the final state, following 
evolution reversal. This overlap is given by 

where Pend(O) is the end-polarized state, and the observ- 
able is the collective magnetization pth = cTz- To low- 
est order, U\jq (see Eq. ([3])) is approximately the in- 



IV. MULTIPLE QUANTUM DYNAMICS: 
BEYOND SPIN-CHAIN APPROXIMATION 



In order to understand the discrepancies observed be- 
tween the analytical model and the experimental results, 
it is necessary to identify the dominant sources of non- 
ideality in the experiment, and assess their respective 
effects on the obscrvablcs under examination. With this 
in mind, we first analyze effects due to limited control, 
such as the higher-order terms in AHT as well as im- 
perfect system initialization. We then investigate the 
intrinsic limitations of the ID NN model to describe 
the real physical system, which contains an ensemble of 
weakly-coupled spin chains with long-range intra-chain 
couplings. In particular, we compare the effects of long- 
range interactions first within a single chain and then 
across different spin chains. Note that while long-range 
couplings have been previously accounted for in a per- 
turbative limit [25j . we resort here to exact numerical 
simulations (Appendix VK\ . while also considering other 
experiment-related sources of errors. 




5 10 15 20 25 30 35 



Number of cycles 

FIG. 3: (Color online) Overlap A for end-polarized initial 
state, Eq. (|10p . Normalized experimental data (stars) and 
numerical results for finite- width pulses (solid line) . The chain 
length in the numerical calculation is A^ = 9, whereas the 
pulse length and inter-pulse delay are w — 0.0075 b ~ Ifis and 
A = 0.0225 b, respectively, which are close to the experimental 
values. The shaded area is the time region explored in the 
MQC experiments. Inset: Absolute value of the difference 
between Jo as calculated from the analytical model or the 
rf-pulsed dipolar Hamiltonian with the initial thermal state 
under different chain lengths, as a function of cycle time. Only 
NN couplings are considered in this case. 
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verse of U^jq. Thus, the overlap A is close to maximal 
for short cycle times. The experimental data is shown 
in Fig. [3l The experimental data were normalized by 
fitting the decay to a normalized Gaussian curve. The 
7r/2 pulse length used was w = 0.93 /iS, whereas the de- 
lay A = 2.9 /iS. In normalized units (the NN coupling 
h 8.3 kHz in practice), this corresponds to Tc ~ 0.72, 
indicating that wc arc well within the regime where the 
contributions of the higher order terms can be neglected. 
Even as A is increased to 7.3 in some of the experi- 
ments, Tc only increases to « 1.64 (in normalized units), 
thus still within the range where higher-order corrections 
are unimportant. 

This is confirmed by numerical simulations, also shown 
in the main panel of Fig. [3l We prepared the end- 
polarized initial state in a matrix form for a system of 
9 spins and evolved the system first forward under the 
DQ sequence with pulses along the x'-axis, then backward 
by using y-pulses. Considering that in practice the DQ 
coupling strength 6^ < 8.3 kHz, and that finite-width 
corrections originate primarily from the second-order av- 
erage Hamiltonian, we expect these corrections to be on 
the order of [hjivY < 6.9 x 10"^. As seen in Fig. O 
the overlap from numerical calculations is flat and close 
to unity, confirming that errors due to finite widths and 
high-orders AHT contributions arc small. Comparison 
with the experimental data suggests that other sources of 
error are likely to be responsible for the long-term decay 
of the overlap [i^. In particular, both rf and static- field 
inhomogcneities can result in imperfect 7r/2-pulscs, lead- 
ing to off-axis and pulse-length systematic errors. The 
latter errors are actually minimized by the 16-pulse se- 
quence thanks to the use of phase alternation [28[ . Fur- 
thermore, transient eff'ects of square pulse always exists 
in pulse-driven experiments. Notice that the MQC data 
of Fig. [2] were measured at relatively short times, t < 0.5 
ms for most of the data. This corresponds to 6 cycles, 
thereby to high values of the overlap. 



2. Initialization 

The basic idea for preparing the end-polarized ini- 
tial state from the thermal state was introduced in p^ . 
Starting from equilibrium, we first rotate the nuclear 
spins into the x-y plane by a tt/2 pulse along a direc- 
tion a. We then allow the system to evolve under the 
dipolar Hamiltonian of Eq. ^ for a time ti {— 30.3 jis 
in the experiment, corresponding to 0.25 in normalized 
units), and finally rotate the spins back to the z-axis by a 
second 7r/2 pulse along the —a direction. During time ti, 
the spins at both ends evolve roughly 1/^2 times slower 
than the internal spins, due to the fact that each of them 
has only one nearest neighbor, while any internal spin has 
two. Let Ua describe evolution under the pulse sequence 
7r/2 |q — ti — 7r/2| -, where in the experiment the pulse 
axis a is phase-cycled through the y- and cc-axes. Given 
that the state at time ti is p{ti) = {1/Na) UaPthUl, 




FIG. 4: (Color online) Left panel: Evolution of the end- and 
central- spin polarizations and fidelity during the preparation 
of the end-polarized state. Chain length = 9, a — x,y, 
and Na = 2. The two time points ti = 0.25 and t'l = 0.42 
are marked with vertical dashed lines. Right panel: MQC 
signal of iV = 19 spin chain with preparation time ti = 0.25. 
Intensities are normalized at every time as Jo -I- 2( J2 + J4) = 
1. DQ Hamiltonian with NN-I-NNN couplings and the ideal 
end-polarized initial state - solid lines; DQ Hamiltonian with 
NN-I-NNN coupUngs and initial state synthesized at ti = 0.25 
- circles for Jo, stars for J2, and squares for J4. 



with Nq, being the number of phase-cycling steps, the 
fidelity of the prepared p(ti) relative to the desired end- 
polarized state is 



Tr[pendP(^l)] 

VTr[pL.]Tr[p2(ii)] 



(11) 



The difference between p{ti) and pend is due to the pres- 
ence of zero quantum coherences which arc generated 
by the dipolar Hamiltonian but are not be removed by 
phase-cycling, with leading contributions from residual 
polarization on spins 2 and N — 1, as well as correlated 
states of the form cr^ {a^_-^a~^-^^ + cTi-if^t+i) [3- The left 
panel of Fig. 3] depicts the time dependence of the fidelity 
and the polarization of the end and the central spins. In- 
terestingly, the time that maximizes fidelity (ii = 0.25) 
does not coincide with the time at which the central-spin 
polarization is zero {t'l = 0.42). It is also worth mention- 
ing that both time points are almost independent of the 
chain length unless iV < 4. 

Starting from the two prepared states, ti = 0.25 and 
t[ = 0.42, respectively, we calculate the MQC of the spin 
chain under the DQ Hamiltonian with NN+NNN cou- 
plings, and compare the results against those obtained for 
the ideal end-polarized initial state PendiO). The evolu- 
tion of MQC for the initial state prepared with t[ = 0.42 
is quite different from that obtained with the intended 
state (data not shown), while the MQC of the initial 
state corresponding to preparation time ti = 0.25 is very 
close, as demonstrated in the right panel of Fig. |4l Note, 
however, that compared to the ideal end-polarized state, 
the experimentally prepared initial state shows slightly 
larger oscillations, especially in Jq. 
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FIG. 5: (Color online) Effect of NNN interactions in the DQ Hamiltonian for the thermal (left column - a and c) and the 
end-polarized initial state (right column - b and d). The length of the spin chain is A'' = 18 (top row - a and b) and 19 (bottom 
row - c and d). Solid (blue and red lines) are Jo and J2 computed from the analytical model, respectively. Circles, stars, and 
squares are Jo, J2, and J4 obtained from exact numerical results. 



B. Non-idealities in isolated single-chain dynamics 

1. Long-range couplings 

As already remarked, appreciable growth of the four- 
quantum coherence signal at long times (Fig. [J) indicates 
the inadequacy of the analytical model, as the dynam- 
ics is no longer confined to a ID system with pure NN 
couplings. Having shown in the previous section that the 
effect of higher-order AHT terms is negligible in the tem- 
poral region we consider, we next turn our attention to 
the influence of long-range couplings. We limit our cal- 
culations to the NNN couplings as they provide the most 
important correction to the analytical NN model. 

Figure [5] shows the MQC signal obtained for spin 
chains of length TV = 18 (top) and 19 (bottom) for the 
thermal state (left) and the end-polarized state (right), 
respectively. Both NN and NNN couplings in the DQ 
Hamiltonian are now exactly accounted for. By way of 
comparison, we also include the predictions from the an- 
alytical model. The following main observations may be 
made: 

(i) NNN couplings produce even-order coherences 
greater than two, the largest contributions in the rele- 
vant time window arising from J4. In general, even order 
coherences up to the number of spins in the chain may 
be expected. This is in contrast with the results based 
on a perturbative approach [2^, which yield MQC only 
up to the sixth order. 

(ii) NNN couplings reduce the amplitude of the oscil- 
lations in Jo and J2. 

(iii) The effect of NNN couplings is amplified at an in- 
stant in time that we call the mirror time, tm (~ 5 in the 
figure) , which is defined in terms of the analytical model 



as the time where Jg shows a second largest oscillation 
for odd A^ or the lowest point for even N. (Note that 
one could also equivalently define tm as the time where 
the second lowest/largest peak of J2 occurs.) This effect 
is prominent in the numerical simulations, where one is 
necessarily constrained to relatively short chains. Quali- 
tatively (see also Appendix B), the spin dynamics has a 
mirror symmetry about the middle spin, which causes the 
signal of specularly located spins to "interfere construc- 
tively" at the mirror time. This picture can also explain 
why the influence of NNN couplings on the dynamics 
of the chosen collective observable is most pronounced 
at this time: even small deviations from the ideal NN 
dynamics arc able to destroy the interferences and can 
produce significant changes in the observed signal. 



2. Chain length distribution 

Since the defects in the FAp sample are non-uniform, 
the spin chain length has a statistical distribution. Ac- 
cording to the so-called random cluster model [2^ , if de- 
fects are distributed randomly in the infinite ID chain 
with a probability (1 — p), the average chain length 
is N = (1 4- p)/(l — p), and the relative fluctuation 
AN/N = -y/2p /(I + p). For a low percentage of de- 
fects, (p ~ 1, N ^ 1), the chain length distribution can 
be reasonably approximated by a uniform distribution 
of chain lengths. Fig. [5] shows the averaged MQC sig- 
nal for an ensemble of chain lengths. Compared to an 
individual spin chain, the ensemble average washes out 
the long-time oscillations but leaves the short-time oscil- 
lations virtually unchanged. Since the concentration of 
defects is low in the actual sample, we expect that this 
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FIG. 6: (Color online) Ensemble average of MQC signal for 
the thermal (left) and the end-polarized initial state (right), 
as resulting from a DQ Hamiltonian with NN and NNN cou- 
plings. Circles, stars, and squares show Jo, J2 and J4, re- 
spectively, with chain length averaged over TV — [17, 21]. For 
comparison, the solid lines are for a distribution of chains with 
average iV = 50 in the random cluster model, evolved under 
the analytical model, whereas the dotted lines are for the zero 
and double quantum intensities for an infinite chain [25j . 



effect will not be important on the time scales explored 
by the current experiments. 



C. Non-idealities due to coupled-chain dynamics 

Due to the 3D nature of the FAp sample, a given spin 
chain of interest ( "central" spin chain henceforth) is cou- 
pled to all other chains in the crystal via the long-range 
dipolar coupling. Since the distance between two spin 
chains in FAp is about three times the distance of two 
NN ^^F spins, the cross-chain couplings have about the 
same strength as the third-neighbor intra-chain coupling 
within a chain. The combined effect is, however, am- 
plified by the presence of several (six) chains surround- 
ing the central spin chain (recall Fig. 1). Furthermore, 
additional weaker contributions arise from more distant 
chains. Overall, the influence of the cross-chain coupling 
can thus be an important source of deviation from the 
analytical model, as we explore next. 

Exactly modeling the influence of all chains on the cen- 
tral one would require us to simulate the quantum dy- 
namics of a macroscopically large number of spins, which 
is clearly beyond reach. To make the problem numeri- 
cally tractable, we thus need to reduce the many-body 
problem to a simpler model that represents as faithfully 
as possible those features of the real dynamics we are 
directly probing. In order to make sensible approxima- 
tions, it is useful to reconsider the origin of the NMR 
signal in more detail. Let Mc be the number of chains 
present in the crystal sample. In the high-temperature 
approximation, the initial density matrix of the whole 
system can be expressed as 



P3d(0) 



Ma 

m—1 



Eqs. (|4])-([5|). Notice that due to its collective nature, the 
experimentally accessible observable can also be written 
as a sum of contributions from distinct chains. For the 
purpose of making contact with a reduced description 
where a single chain is singled out as a reference in the 
presence of other coupled chains, it is useful to view the 
total signal in the k-th measurement as originating 
from two terms, S^{t) = S'^ ^^^ 



S,.leak' with 



Antra (0 



z.leak 



it) 



E 



Tr 



(12a) 



(12b) 



where m indexes the chains and pm (0) is either the ther- 
mal equilibrium state or the end-polarized state, as in 



where pf„(t) = Uk{t)pm{0)Ul{t) (cf. Eq. 0). These 
two terms reflect two different mechanisms by which the 
presence of cross-chain interactions can induce deviations 
of the experimental signal from that of an isolated chain. 

The term in Eq. (|12ap . which we refer to as the intra- 
chain signal S^^^^^^{t), describes the signal obtained 
when the initial state and observable belong to the same 
chain: all other chains, which may initially be taken to 
be in the maximally mixed state, influence the reference 
chain in a "mean-field sense," to the extent they modify 
Pm{t)- Were all the chains identical, the resulting sig- 
nal would simply be 5'^_j„tra(0 ~ A4Tr[S'zp''(t)], that 
is, an A/c-fold signal from a single chain coupled to the 
"environment chains". Thus, 5'^i„ira(0 ^^^Y be well de- 
scribed within a "chain-plus-environment" model, where 
a single spin chain is coupled to a larger spin environ- 
ment, and the measured NMR signal is determined com- 
pletely by the reduced density matrix of the central chain 
- upon tracing over all environment spins, as in a stan- 
dard formulation of the central-system-plus-bath prob- 
lem [4l,[5Sp. 

While the intra-chain term describes a deviation from 
the ideal single-chain behavior that is not fundamen- 
tally different from, say, deviations induced by long-range 
couplings as analyzed in Sec. IIVB 11 the leakage signal 
'^z leaki^) in Eq- (|12bp introduces a qualitatively different 
effect: that is, the possibility that some of the polariza- 
tion initially located on the mth chain is transferred to 
the m'th chains, and read out there. Since, from the 
point of view of the central spin chain, signal compo- 
nents would be 'lost' to the environment, a significant 
contribution j^^j, (t) would clearly indicate the inade- 
quacy of a systcm-plus-cnvironment picture at capturing 
the complexity of the underlying 3D strongly-correlated 
dynamics. 

Even assuming that the consistency of a central chain- 
plus-environment treatment may be justified a posteriori 
by the smallness of the leakage signal for the evolution 
times of interest, modeling a realistic environment re- 
mains non-trivial because the actual crystal consists of a 
large number of quantum spin chains, evolving according 
to a highly complex, non-Markovian dynamics. In line 
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with standard statistical approaches (including NMR re- 
laxation theories) [28l . [s^ . we can however reasonably 
argue that the main observed features should be robust 
with respect to the details of the environment descrip- 
tion, as long as the relevant energy scales are correctly 
reproduced. In what follows, we will exemplify these con- 
siderations by separately investigating two models for de- 
scribing how the coupling between different chains in FAp 
modifies the MQC dynamics of the central spin chain. In 
Sec. lIVCTl a system of two coupled chains is investigated 
as a numerically accessible testbed where the 'environ- 
ment chain' qualitatively retains the spatial structure of 
the FAp crystal. Physically, the latter feature is expected 
to be important (possibly essential) to properly represent 
the deviation induced in the idealized central-chain dy- 
namics by the NN chains. In Sec. IIV (J2l a structureless 
spin environment model is considered instead, whereby 
the central chain couples to randomly placed spins. Phys- 
ically, such a picture may be especially adequate to ac- 
count for the net influence of distant chains. Computa- 
tional constraints limit the size of the accessible model 
environment in both cases. 

In spite of the above differences, it is important to re- 
alize that essentially the same type of simulations will be 
employed and the same main physics will be explored in 
both cases. In particular, the processes leading to devia- 
tions from the analytical model are primarily associated 
with the increased dimensionality of the Hilbert space 
and correlations between different chains in the sample. 
While no explicitly non-unitary evolution is present ci- 
ther in experiment or simulation, and the total system 
remains coherent at all times, a damping of the low- 
order MQC oscillations still emerges: as time progresses, 
a larger part of the Hilbert space is populated, and co- 
herences of higher-order, which involve spins of different 
chains, build up at the expenses of low-order coherences. 
Relative to the observables that can be directly probed, 
the latter simply appear to unrecoverably decay. 



1. Effect of a structured spin environment 

The correlated dynamics in nearby chains may be in- 
vestigated by lumping together the contributions of the 
six nearest surrounding chains and treating them as a 
single chain, which couples coherently to the central spin 
chain, according to the DQ Hamiltonian. We take both 
chains to have length N and start in the initial state of 
interest (either thermal or end-polarized). Since we are 
restricted to numerically calculate MQC for a system of 
up to 25 spins, A'^ < 12 in practice. Upon summation 
[Mc = 2 in Eq. p2a| - p2b|) ]. an upper bound to the NN 
cross-chain coupling strength is given by 




-0.1488. 



(13) 



This approximation corresponds to neglecting correla- 
tions between spins from three or more different chains, 



FIG. 7: (Color online) Effect of cross-chain couplings for the 
thermal (top) and the end-polarized state (bottom). The 
length of each spin chain is N = 11. Dotted lines represent 
the analytical isolated-chain NN prediction. Both NN and 
NNN couplings are considered in the exact simulation results 
for a single chain (dashed lines) as well as for two chains (solid 
lines) coupled according to Eq. (|13|l . 



which arise from higher-order cross-chain couplings in 
leaki^) ^ ^'^^ instance, the three-chain coupling is pro- 
portional to (bx/b)^- If such couplings are treated per- 
turbatively, one may expect their effect to be negligi- 
ble over the time scale of the experiment, as opposed to 
two-chain interactions which directly compete in strength 
with intra-chain NNN couplings. As discussed above, 
however, these two contributions may have very differ- 
ent physical implications, as the cross-chain coupling ef- 
fect can genuinely increase the underlying Hilbert space, 
whereas NNN couplings can only increase the portion of 
the single spin-chain Hilbert space that is explored dur- 
ing the dynamics. 

Exact calculation of the total signal (t) reveals that 
the contribution of cross-chain transfer due to i^^f. (t) 
remains small (below a few percents) over relatively short 
time scales (up to 5 in normalized units). As shown in 
Fig. [71 cross-chain couplings modeled in this way also 
damp the MQC oscillations at long times, similar to the 
effect of intra-chain NNN interations. Notice that at the 
mirror time, as observed in simulations with finite N, 
the effects of the cross-chain couplings are also amplified, 
further reducing the peak amplitude [53| . 



2. Effect of a structureless spin environment 

According to our earlier discussion, another way to an- 
alyze the influence of cross-chain coupling that empha- 
sizes the influence of far away chains is to consider an 
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FIG. 8: (Color online) Jo{t), Jiit), and J4,{t) for a chain of 
TV = 11 spins with thermal (left) and end-polarized (right) 
initial state. The chain is coupled to an environment of nine 
spins with random dipolar couplings. The lines are the ana- 
lytical results for Jo(f), J2{t). Note that only NN intra-chain 
couplings are included, however the times scales explored here 
are significantly longer than in any of the previous figures. 



effectively structureless quantum spin environment [5J|. 
In particular, the simplest choice is provided by a sys- 
tem consisting of spins arbitrarily (randomly) scattered 
in space. While of course nothing is arbitrary in the dy- 
namics of the real FAp system, such a randomized model 
may just be viewed as a computationally accessible ap- 
proximation of the complex dynamics under investiga- 
tion. 

Specifically, we reproduce the main features and the 
characteristic energy scales of the FAp sample driven by 
the 16-pulse sequence by assuming that the x, y, z co- 
ordinates of each of the 9 environment spins are drawn 
uniformly from [—1,1]. The = 11 spins of the cen- 
tral chain are placed equidistantly on the z-axis, with 
their z-coordinates also confined between —1 and 1. The 
minimum distance between any pair of spins (whether 
environment or chain spins) is restricted to exceed 0.1 to 
prevent spins from being too close to each other. The 
central chain Hamiltonian of the form ([1]) is truncated at 
either the NN or NNN level (Fig. [8] and [9] respectively) . 
All the dipolar coupling between the environment spins, 
and from the environment spins to the central spin chain 
are taken into account, as in Eq. (IT|), with the coupling 
constants hjn calculated from the spins coordinates. How- 
ever, in order to have correct energy scales, all chain- 
environment coupling constants are rescaled to produce 
the correct value of the dispersion Tr where Hcb 

is the chain-environment interaction Hamiltonian (see 
also Appendix D). This ensures that the couplings be- 
tween the spins of different chains in FAp are « 40 times 
smaller than the couplings between the spins in the same 
chain. In a similar way, all couplings inside the envi- 
ronment are rescaled to produce a correct value for the 
Hamiltonian norms per spin, Tr [iy^"] ~ \^'b\ ' where 
_ff(7, Hb are the chain and environment Hamiltonians, 
respectively. 

We perform simulations of the total system treating it 
as a closed quantum system with unitary dynamics (see 
Appendix [X] for details). We simulate the evolution un- 
der the experimental DQ Hamiltonian, generated by the 




FIG. 9: (Color online) Comparison between the effect of cross- 
chain couplings starting from the thermal state, as resulting 
from a structured two-chain system (solid lines) or from a cou- 
pling to a randomized spin (circles). In all cases, the length 
of the spin chain(s) is = 11 and NNN intra-chain couplings 
are included. Dotted lines represent the analytical predic- 
tion. Numerical results for Jo(t) (blue), J^it) (red), and J^it) 
(black) MQCs are included. 



16-pulse sequence. Each sequence with 16 pulses along 
the X-axis (preparation) and, afterward, the 16-pulse se- 
quence with pulses along the y-axis is repeated five times 
(mixing). The pulses are ideal (5- like, with varying inter- 
pulse separation, and the total time is varied from zero to 
18.75 (in normalized units). Note that since the environ- 
ment is homonuclear, it is affected by the pulses in the 
same way as the central chain. At the end of each proto- 
col, the total NMR signal (t) is calculated by either (i) 
summing only the z-projections of the spins in the central 
chain, tracing out the environment spins (thus obtaining 
only the intra-chain contribution of Eq. (|12a|) : or (ii) 
summing the z-projections of both the environment and 
the chain spins (thus also taking into account the leakage 
terms in Eq. (jl2bp ). Comparison between the results (i) 
and (ii) shows that the leakage terms are small, on the 
order of about 1%. As in the two-chain model, we thus 
confirm a posteriori the validity of the underlying weak- 
coupling assumption between the central system and the 
rest. 

Numerical results starting from the thermal and the 
end-polarized state are given in Fig. [S] for a single real- 
ization of such a random dipolar ly-coupled environment, 
corresponding to a fixed (arbitrary) geometry of the spin 
lattice. While different realizations give very close results 
(data not shown), averaging over several realizations is 
impractical [s^. Fig. [8] also includes a comparison of the 
simulation results for Jo. 2(0 with the analytical model. 
The interaction with the environment leads to a signifi- 
cant damping of the oscillations of Jo{t) and J2(i), and 
to an overall decay of these coherences. Interestingly, the 
decay of both Jo (t) and J2 (t) for the end-polarized initial 
state is slower than for the thermal initial state. Likewise, 
it is also worth noticing that the decay of the oscillations 
in J2 is roughly a factor of two slower than the decay of 
Jo. This difference may be attributed to the fact that 
the random dipolarly-couplcd environment is not fully 
structureless, as it possesses non-trivial integrals of mo- 
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tion (for instance the total magnetization of the central 
chain and the environment). The internal structure of 
the environment appears to strongly affect the dynamics 
of J2 (<) . This different behavior of the two MQC intensi- 
ties is also present in the two coupled-chain simulations 
of Sec. IIV C 1\ which are directly contrasted to the ran- 
dom spin-environment simulation results in Fig. [S] We 
further expand on these considerations by examining a 
chaotic spin bath model in Appendix [Cl 



V. DISCUSSION AND CONCLUSION 

We have investigated in detail the MQC dynamics of 
a quasi-lD spin chain in a fluorapatite crystal, both ex- 
perimentally and numerically. By comparing exact sim- 
ulation results with analytical solutions for the ideal DQ 
Hamiltonian with NN couplings, we have characterized 
the region of validity of this simple, single-chain NN 
model. For the initial states and observables of inter- 
est, we have found that for evolution times up to 0.5 ms 
(corresponding to about 5 times the inverse NN coupling 
strength) the system is experimentally indistinguishable 
from the single-chain, NN model. Simulations including 
long-range couplings within a single chain and across dif- 
ferent chains reproduce well the experimental findings. 

Beyond this time, the evolution deviates from the an- 
alytical model, although the deviations of the selected 
observables (the MQC) remain small. In principle, the 
experimental implementation of the DQ Hamiltonian us- 
ing a simulation approach based on AHT is not a prob- 
lem, at the evolution times considered. In addition, the 
dynamics of the experimentally created end-polarized ini- 
tial state are seen to remain quite close to the dynamics 
of an ideal end-polarized state, as desired. 

From simulations we observed that all the different 
types of long-range couplings analyzed lead to a quali- 
tatively similar damping of the oscillations in the MQC 
signal and a relatively slow growth of the higher order 
coherences (in particular the 4-quantum coherence). In 
fact, a similar effect is also observed for a single chain 
coupled to a dipolar spin environment. 

The similar behavior observed when introducing longer 
range couplings in a ID chain and cross-chain couplings 
seems to indicate that although in the second case there 
are more pathways available for the propagation of multi- 
spin correlations, this effect cannot be observed in the 
MQC evolution. While it could be tempting to infer that 
the microscopic mechanisms leading to the observed be- 
havior are to some extent similar in each case, it is also 
essential to acknowledge that the experimentally accessi- 
ble, collective magnetization observable provides a highly 
coarse-grained visualization of the overall dynamics. 

From a many-body physics standpoint, a deeper un- 
derstanding of the influence of the structure of the 
longer-range dipolar couplings ( "internal environment" ) 
on MQC dynamics, in particular of the potentially higher 
level of sensitivity found for higher-coherence orders, is 



certainly very desirable. 

Lastly, from a quantum communication perspective, 
our work calls attention to the added challenges that 
transport protocols need to face in the presence of lim- 
itations in available control, initialization, and readout 
capabilities, as well as long-range interactions and/or 
unwanted interactions with uncontrolled degrees of free- 
dom. Our study points out that for the realization of 
precise transport, simply isolating a ID system is not 
enough, as the deviation from an ideal NN model in a 
ID chain caused by long-range couplings is as impor- 
tant as cross-chain couplings. Since a number of these 
issues are shared by all practical device technologies to a 
greater or lesser extent, it is our hope that our analysis 
will prompt further theoretical investigations of commu- 
nication protocols under realistic operational and physi- 
cal constraints. 
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APPENDIX A: NUMERICAL METHODS 

We calibrate our numerical procedure by reproducing 
the results from the analytical model [l^, [l3| . For suffi- 
ciently short spin chains, < 11, we propagate exactly 
the density matrix of the system. That is, given an ini- 
tial mixed state (either the thermal or the end-polarized 
state), we prepare the initial density matrix, evolve the 
system, obtain the density matrix at time t, and calcu- 
late the MQC signal according to Eq. ([7]). For longer spin 
chains, this approach becomes very inefficient due to the 
extremely high usage of computer memory (on the order 
of 4^ for a chain with N spins). Instead, we employ a 
wave-function-based simulation method, whose memory 
usage scales as 2^. 

To implement the wave-function simulation, we decom- 
pose the initial mixed density matrix of the system into 
a sum of (for the thermal state) or 2 (for the end- 
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polarized state) individual density matrices, and then ap- 
proximate the jth density matrix with a product state of 
a known state of spin j and a pure random state of the 
remaining spins [56l. |57|. That is, we let 

'^]-l{h-\i^){i,\)®h){r,l 

where \rj) = X]i=i is a linear combination of ba- 

sis states of all spins except the jth spin, and is a 

random complex number obeying X]i=i ~ 1- Such 

a superposition is an exponentially accurate representa- 
tion of the maximally mixed state, and in our simula- 
tions creates errors of about 0.5%. After preparing the 
initial wave- function, we propagate the system according 
to the Schrodinger equation, adopting an efficient algo- 
rithm based on Chebyshev polynomial expansion of the 
evolution operator |58| . 

In the calculation of the MQC signal for the system- 
plus-environment, an alternative way to prepare the ini- 
tial state is used, by realizing that the initial density 
matrix may be expressed in terms of spin operators as 

follows. Let \R) = "^iK) ^"-^ ^ random wave-function 

of N spins, and \R') = X^jLi Then we may sim- 

ply write p(0) = \R'){R\. The propagation of these two 
wave- functions is then implemented based on the meth- 
ods mentioned above. 
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FIG. 10; (Color online) (a) Jo for the end-polarized initial 
state of different length spin chains (DQ Hamiltonian with 
NN-fNNN couplings). The dashed line is a visual guide to 
the mirror time, (b) Jo for the thermal initial state and (c) 
its partition for different locally polarized initial state erf (i — 
1, 2, ■ ■ ■ , 9). Due to the mirror symmetry of the chain, only 
half of the initial states are presented in panel (c). Each curve 
is shifted upward 0.5 for a better view. The vertical dashed 
line specifies the position of mirror time tm- For (b) and (c), 
the length of the spin chain is A'^ = 17. 



APPENDIX B: MIRROR TIME 

Besides the peak (dip) of the MQC signal and the am- 
plification of the NNN coupling effects at the mirror time 
t„i, the following features may be interesting for spin 
transport in short spin chains: 

(i) The mirror time increases linearly with the length 
of the spin chain A^, as shown in Fig. flOT a). 

(ii) For the same length spin chain, different locally po- 
larized initial states have the same mirror time (of course, 
the thermal state, which may be seen as a mixture of dif- 
ferent locally polarized initial states, also exhibits the 
same mirror time); See Fig. llOf b)-(c). 

(iii) The NNN couplings shift the mirror time slightly. 

These peculiar properties demand a better understand- 
ing of the physical meaning of the mirror time. In a 
picture of spin polarization transport along a chain p7| . 
starting from the end-polarized state where the polar- 
ization is pinned to spins 1 and A^, the polarization is 
transported to the central spin (A — l)/2 at the mir- 
ror time tm (we assume A^ is odd for simplicity). As 
mentioned in the main text, the spin dynamics exhibits 
a mirror symmetry about the central spin, and thus in- 
terferes constructively at i,„. For other pairs of locally 
polarized initial states, for instance (t| and af^^^_j, the 
spin polarization also interferes constructively at the mir- 
ror time. The independence of on j guarantees that 
the thermal state shows the same properties at t„j as the 
end-polarized state. 



APPENDIX C: CHAOTIC BATH MODEL 

Since in simulations we cannot exactly reproduce the 
many-body dynamics occurring in the FAp crystal, ap- 
proximations are necessary at a number of levels. In 
representing the dynamics in terms of a single chain cou- 
pled to a bath the random dipolarly-coupled environment 
model used in the main text (Sec. lIVC2|) imposes a struc- 
ture on the environment that is motivated by the physical 
system itself. From an open-system perspective, how- 
ever, it may be interesting to explore alternative models 
for the bath, in order to have a sense of which details are 
important for the system's dynamics and which are not. 
Although these alternative bath models need not have 
an immediate relevance to the experimental system, they 
may provide additional physical insight on the action of 
a spin bath in the FAp crystal. In this venue, it is useful 
to observe that quantum systems possessing a very com- 
plex behavior often exhibit similar features, and relevant 
aspects of their dynamics may be captured by quantum 
chaotic models, see e.g. [s^. Following this approach, 
we emulate the bath's internal dynamics using a chaotic 
spin-glass s/iard Hamiltonian [60l. [6l|. As a main feature, 
the chaotic bath model assumes that no integrals of mo- 
tion exist for the bath other than the energy. This differs 
from the dipolarly-coupled environment model (and the 
real FAp sample), where the environment chains are sim- 
ilar to the central chain and, in the absence of pulses, the 
total magnetization of the central chain and the bath is 
conserved. 
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FIG. 11: (Color online) Jo(t), J2{t), and J4{t) for a chain 
with thermal initial state (left), and for the chain with end- 
polarized initial state (right). The chain is coupled to the bath 
of nine spins with chaotic glass-shard internal dynamics de- 
scribed by Eq. (jCip . Dash-dotted lines: Simulations with the 
dipolar ly-coupled spin environment described in Sec. IIV C 21 

Specifically, in our case we choose the chain-bath cou- 
pling to mimic the arrangement of FAp samples: each 
chain spin is coupled to six bath spins, the coupling has 
a homonuclear secular dipolar form, similar to Eq. p]). 
and the coupling constants bj£ for each pair of a chain 
and a bath spin are drawn uniformly from the interval 
[-V3 • 0.025, ■ 0.025]. This ensures that the rms cou- 
pling between one bath spin and one chain spin is equal 
to the experimental value bx/b w 0.025, see Eq. ([T5|. 
Nine bath spins are located on a 3 x 3 square lattice, 
with a Hamiltonian 

(fe,i) k k 

where the summation in the first term is over NN pairs. 



The random couplings Tki and the local magnetic fields 
h^'^ are drawn uniformly from the intervals [— Eg, Eq] and 
[~ho, ho], respectively, with the values of Eq and ho ad- 
justed to ensure: (i) chaotic regime, and (ii) correct char- 
acteristic energies for the spin dynamics inside the bath 
To achieve the latter, note that for a EAp chain with NN 
couphngs only, and TV > 1 spins, TtH^ = (6/16) N Trl, 
so that the rms energy per spin is 6/16. Correspondingly, 
the values of Eq and ho were adjusted to give approxi- 
mately the same rms energy per spin. 



The results of the simulations for the thermal initial 
state and for the end-polarized initial state are given in 
Fig. [TlJ It is clearly seen that the interaction with the 
bath leads to significant damping of the oscillations of 
Jo{t) and J2{t), and to an overall decay of these coher- 
ences, although the mirror time remains clearly visible. 
Interestingly, as also noted in the text, the decay of the 
zero- and second-order coherences Jo{t) and J2{t) for the 
end-polarized state is slower than for thermal state. To 
further appreciate this, we compare the dynamics of Jo(t) 
and J2{t) for the two bath models we examined in Eig.fTT] 
The Jo{t) signals for both bath models stay close to each 
other, while exhibiting significant damping of oscillations 
and overall decay in comparison with the analytical re- 
sults for the isolated chain. In contrast, J2{t) for the ran- 
dom dipolarly-couplcd environment stays rather close to 
the analytical prediction for the isolated chain, whereas 
J2 (i) for the chaotic bath decays in the same way as Jq (t) 
does. This suggests that the presence of extra integrals of 
motion does not significantly affect the dynamics of Ja{t), 
whereas higher-order MQCs might more sensitively de- 
pend upon details of the open-system dynamics. 
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